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Abstract 

The MiniBooNE Collaboration observes unexplained electron-like events in the reconstructed 
neutrino energy range from 200 to 475 MeV. With 6.46 x lO^'^ protons on target, 544 electron- 
like events are observed in this energy range, compared to an expectation of 415.2 it 43.4 events, 
corresponding to an excess of 128.8 it 20.4 zb 38.3 events. The shape of the excess in several 
kinematic variables is consistent with being due to either Ve and Ue charged-current scattering or 
to neutral-current scattering with a photon in the final state. No significant excess of events is 
observed in the reconstructed neutrino energy range from 475 to 1250 MeV, where 408 events are 
observed compared to an expectation of 385.9 it 35.7 events. 
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In a previous Letter [l| , the MiniBooNE collaboration reported initial results on a search 
for — > Ue oscillations. The search was motivated by the LSND observation ^ of an 
excess of Ue events in a beam that implied larger values of Am^ than any of the currently 
confirmed oscillation measurements. The MiniBooNE result showed no evidence of an excess 
of electron-like events for neutrino energies above 475 MeV. However, a sizeable excess 
of electron-like events was observed from 300-475 MeV. This Letter reports on a more 



detailed investigation of the low-energy electron-like events [31]. Published explanations for 
the low-energy excess range from anomaly jmediated neutrino-photon coup ling J] to neutrino 
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oscillations involving sterile neutrinos [5|, 16|, |7|, 18|, |9|] to Lorentz violation [10|] . In the course 
of this investigation, many improvements have been made to the data analysis, and the data 
sample has increased from 5.58 x 10^° protons on target (POT) to 6.46 x 10^*^ POT. The 
excess of electron-like events persists after these improvements and has been studied as a 
function of several kinematic variables. 

MiniBooNE uses the Fermilab Booster neutrino beam, which is generated from 8-GeV 
kinetic energy protons incident on a beryllium production target. Neutrinos are produced 
in a 50 m long decay pipe by the in-fiight decay of pions and kaons and a small fraction of 
the sub seq uent muons. The center of the MiniBooNE detector is 541 m from the production 
target ll|. The neutrino target and detector medium is mineral oil in which relativistic 
particles create both Cherenkov and scintillation light. The different properties of these 
sources of light readily allow particle identification; however, the detector cannot distinguish 
between electrons and photons. 

The Booster neutrino beam flux at the detector is modeled using a GEANT4-based 
simulation Il2| of the beamline. Pion and kaon production in the target is parametrized 



by a global fit to proton-beryllium particle production data 



The z/o energy 



spectrum peaks at ~ 600 MeV, has a mean energy of ~ 800 MeV, and extends to ~ 3000 
MeV H. 

The specific changes to the analysis of the low-energy events since the initial paper |l| 
are discussed in some detail in the following text. 



The v3 NUANCE 



17| event generator is used to simulate neutrino interactions in min- 



eral oil. The constraint on neutral-current (NC) 7r° production from MiniBooNE data was 



expanded to finer momentum bins 



Also, a direct measurement of low energy NC coher- 



ent vr production was implemented to improve the modeling of n events in the most 
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forward direction. In addition, there is a more accurate treatment of the ratio of 7 to 7r° 
decay of A in nuclei. To avoid uncertainties in neutrino flux and NC cross sections, the 
number of A radiative decays is determined from the number of measured NC 7r° events. 

n 

Final state particles from the initial neutrino interaction [17|, their decays, and possi- 
ble strong and electromagnetic re-interactions in the detector medium are modeled using 
a GEANT3-based 19| simulation, with strong interactions simulated using GCALOR 20|. 



Since the previous Letter Ij, a number of processes, missing from the strong interaction 
model, have been added that could create electron-like backgrounds: photonuclear interac- 
tions on carbon, radiative vr" capture, radiative decay of A resonances produced in pion- 
carbon interactions, and vr^-C (strong) elastic scattering. Radiative capture and A N'j 
decay produce single photons that MiniBooNE cannot distinguish from electrons. Photonu- 
clear interactions can cause a photon from a vr" to be missed, leaving a single photon. Elastic 
scattering of charged pions can cause Cherenkov rings to appear more electron-like. Of these, 
only photonuclear interactions contribute significantly to the electron-like background with 
apparent neutrino energy > 200 MeV. The well-measured photonuclear cross section on 
carbon is used to simulate final states from excitation of the giant dipole resonance and A 
production above and below the pion threshold. The addition of photonuclear absorption 
increases the estimated background from NC vr^ scattering by ~ 30% in the energy range 
200 < E^^ < 475 MeV. E^^ is the reconstructed incident neutrino energy and is deter- 
mined from the reconstructed lepton energy and angle with respect to the known neutrino 
direction, assuming charged-current quasi-elastic (CCQE) scattering. 

One of the larger Uc backgrounds at low energy results from neutrino interactions in the 
tank wall and concrete vault and dirt surrounding the detector. These events originating 
outside the detector are uniquely characterized by low visible energy (Ems), large radius, 
and a direction that points into the detector; therefore, their contribution can be measured 
from MiniBooNE data. An improved estimate of this background using reconstructed event 
position and direction information reduces the normalization of such backgrounds by 30%. 
In addition, a new selection criterion based on energy and topology rejects 83% of these 
events, while discarding only 21% of signal events in the 200 < E^^ < 475 MeV energy 
range. 

Numerous improvements have been incorporated in the systematic error determination 
associated with the neutrino flux, detector, and neutrino cross section modeling. In estimat- 
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ing neutrino flux uncertainties, the propagation of production errors has been upgraded 
to remove unnecessary model dependence. This results in a decrease in the overall tt"^ pro- 
duction uncertainty from ~ 16% to ~ 7% jl6|, which better reflects the uncertainties in the 
underlying HARP measurement of 7r+ production on Be jl4|. In the detector simulation, a 
comprehensive set of flnal state variations has been evaluated to conservatively encompass 
the uncertainty in the aforementioned list of added hadronic processes. These uncertainties 
contribute an additional 1% uncertainty in the low energy MiniBooNE oscillation analysis. 
In the neutrino cross section model, the estimation of the A radiative decays uncertainty 
has increased from 9% to 12%. Also, measurements of the rates of coherently and resonantly 
produced vr" events |18|] has enabled some reduction in these errors. 

The reconstruction and selection of electron-like events is identical to the initial analysis 
[l| with the addition of the cut to reject events produced outside the detector described 
earlier. Events are reconstructed under four hypotheses: a single electron-like Cherenkov 
ring, a single muon-like ring, two photon-like rings with unconstrained kinematics, and two 
photon-like rings with an invariant mass M^^ = 771^^0. To select z/g-candidate events, an 
initial selection is flrst applied followed by particle identification cuts. 

Four different analyses are performed on the data. 

• Original Analysis: original analysis [l| with the original data set of 5.58 x 10^° POT. 

• Revised Analysis: the Original Analysis with the updated background and uncertainty 
estimates described in this paper. 

• Extended Analysis: the Revised Analysis but with the extended data set of 6.46 x 10^° 
POT. 

• Final Analysis: the Extended Analysis but including the new external event cut. 

Table [I shows the expected number of events with E'^^ between 200 — 300 MeV, 300 — 475 
MeV, and 475 — 1250 MeV after the complete event selection of the Final Analysis. The 
background estimates include antineutrino events, representing < 2% of the total. The total 
expected backgrounds for the three energy regions are 186.8 ± 26.0 events, 228.3 ± 24.5 
events, and 385.9 ± 35.7 events, respectively. 

A total of 1069 events pass the complete event selection of the Final Analysis with 
E^^ > 200 MeV. The numbers of data, background, and excess events for different E^^ 



TABLE I: The expected number of events in the 200 < E^^ < 300 MeV, 300 < E^^ < 475 MeV, 
and 475 < E^^ < 1250 MeV energy ranges from all of the backgrounds after the complete event 
selection of the Final Analysis. 
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186.8 ± 26.0 
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ranges are shown in Table [TTl together with the significance of the excesses for the four 
analyses. The uncertainties include both statistical and systematic errors. While there is 
no significant event excess for E^^ > 475 MeV, a sizeable excess is observed for E^^ < 475 
MeV. For the Final Analysis, an excess of 128.8 ± 20.4 ± 38.3 events (3.0cr) is observed for 
200 < E^^ < 475 MeV. 

Figure [T] shows the E^^ distribution for data (points with statistical errors) and back- 
grounds (histogram with systematic errors) for the Final Analysis, and Fig.[2]shows the event 
excess as a function of E^^. Also shown in the figure, for comparison, are expectations from 
the best oscillation fit and from neutrino oscillation parameters in the LSND allowed region 
fl, which are ruled out at 95% CL if the data are fit with E^^ > 475 MeV The error 
bars include both statistical and systematic errors. The best ocillation fit for E^^ > 200 
MeV corresponds to Am^ = 3.14 eV^ and sin^ 26^ = 0.0017 and has a x'^/DF = 18.3/17. 
The null fit has a x'^/DF = 22.0/19. For > 475 MeV, the best fit is consistent with 
the initial result of no oscillations |l| . As shown in Fig. [3] for E^^ > 200 MeV, the event 
excess occurs for E^s < 400 MeV. 

Figs, m and [5] show the event excess as functions of reconstructed and cos(6') for 
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TABLE II: The number of data, background, and excess events for different E^^ ranges, together 
with the significance of the excesses. The different analyses are described in the text. 



Event Sample Original Analysis [Ij Revised Analysis Extended Analysis 


Final Analysis 


200 - 300 MeV 

Data 375 368 427 
Rarkp-rmind 98*^ + 37 339 4 + '^8 Q 38fi D + 44 3 

Excess (Significance) 92 ± 37 {2.5a) 35.6 ± 38.9 (0.9cr) 41.0 ± 44.3 {0.9a) 


232 

45.2 ± 26.0 (1.70-) 


300 - 475 MeV 

Data 369 364 428 

Rarlrcrrniinrl 97*^ + 9fi 9S9 Q + 98 '^'^0 D + "^1 R 

Excess (Significance) 96 ± 26 (3.7a) 81.1 ± 28.3 (2.9cr) 98.0 ± 31.8 (3.1a) 


312 

83.7 ±24.5 (3.4cr) 


200 - 475 MeV 

Data 744 732 855 
Background 556 ± 54 615.3 ± 58.0 716.1 ±66.2 
Excess (Significance) 188 ± 54 (3.5cr) 116.7 ± 58.0 (2.0cr) 138.9 ± 66.2 (2.1cr) 


544 
415.2 ±43.4 
128.8 ±43.4 {3.0a) 


475 - 1250 MeV 

Data 380 369 431 
Background 358 ± 40 356.0 ± 33.3 412.7 ±37.6 
Excess (Significance) 22 ± 40 (0.6c7) 13.0 ± 33.3 (0.4o-) 18.3 ± 37.6 (0.5c7) 


408 
385.9 ±35.7 
22.1 ±35.7 (O.60-) 


300 < E^^ < 475 MeV, the energy region with the most significant excess. is determined 
from the energy and angle of the outgoing lepton, assuming CCQE scattering, and 6 is the 
angle between the incident neutrino and outgoing lepton. Also shown in the figures are the 



expected shapes from the NC vr'' and A — > N'y reactions, which are representative of photon 
events produced by NC scattering, and from z/gC e~X and UeC e^X CC scattering. 
The different reactions all assume the same u^^ energy spectrum. As shown in Table IIIIl the 
values from comparisons of the event excess to the expected shapes are acceptable for all 
of the processes. Also shown in the table is the factor increase necessary for each process to 
explain the low-energy excess. In each case, the estimated background would have to more 
than double (increase by > 5a) to explain the excess. 
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FIG. 1: The E^^ distribution for data (points with statistical errors) and backgrounds (histogram 
with systematic errors). 
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FIG. 2: The event excess as a function of E^^ . Also shown are the expectations from the best 
oscillation fit and from neutrino oscillation parameters in the LSND allowed region [3]. The error 
bars include both statistical and systematic errors. 
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FIG. 3: The event excess as a function of Eyis for E^^ > 200 MeV. Also shown are the expectations 
:'rom the best oscillation fit and from neutrino oscillation parameters in the LSND allowed region 
^. The error bars include both statistical and systematic errors. 
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FIG. 4: The event excess as a function of for 300 < E'j}^ < 475 MeV. Also shown are the 
expected shapes from the NC vr*^ and A Nj reactions, which are representative of photon events 
produced by NC scattering, and from CC u^C — > e~X and v'eC e^X scattering. The error bars 
include both data statistical and shape-only systematic errors. 



S 50 

« 40 
0) 

u 

uj 30 
20 



data - expected background 
7i" background shape 
A background shape 
Vj, signal shape 
signal shape 




FIG. 5: The event excess as a function of cos(^) for 300 < E^^ < 475 MeV. The legend is the 
same as Fig. HI 

In summary, MiniBooNE observes an unexplained excess of 128.8 ± 20.4 ± 38.3 electron- 
like events in the energy region 200 < E^^ < 475 MeV. These events are consistent with 
being either electron events produced by CC scattering {ugC e~X or v^C — e'^X) or 
photon events produced by NC scattering {vC w-fX). Upcoming MiniBooNE results 



with the Booster antineutrino beam and with the NuMI neutrino beam 



211 ] should help 



distinguish these two possibilities and shed further light on the low-energy region. 
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TABLE III: The values from comparisons of the event excess and cos{9) distributions for 
300 < E^^ < 475 MeV to the expected shapes from various NC and CC reactions. Also shown is 
the factor increase necessary for the estimated background for each process to explain the low-energy 

excess. 

Process x^(cos6')/9 DF X^(Q^)/6 DF Factor Increase 

13.46 2.18 2.0 

16.85 4.46 2.7 

14.58 8.72 2.4 

10.11 2.44 65.4 



the A -H- A^7 background. 



NC 
A ^ iV7 
VeC e-X 
UeC e+X 
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